Highly Entangled Ground States in Tripartite Qubit Systems 



Beat Rothlisberger, Jorg Lchmann, D. S. Saraga, Philipp Trabcr, and Daniel Loss 
Department of Physics, University of Basel, Klingelbergstrasse 82, CH-4056 Basel, Switzerland 

(Dated: February 7, 2008) 

We investigate the creation of highly entangled ground states in a system of three exchange- 
coupled qubits arranged in a ring geometry. Suitable magnetic field configurations yielding approx- 
imate GHZ and exact W ground states are identified. The entanglement in the system is studied at 
finite temperature in terms of the mixed-state tangle r. By generalizing a conjugate gradient opti- 
mization algorithm originally developed to evaluate the entanglement of formation, we demonstrate 
that r can be calculated efficiently and with high precision. We identify the parameter regime for 
which the equilibrium entanglement of the tripartite system reaches its maximum. 
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Entangled quantum systems have been the focus of 
numerous theoretical and experimental investigations 
U B El- In particular, entanglement has been identi- 
fied as the primary resource for quantum computation 
and communication 0]. Compared to the case of a bi- 
partite system, multipartite entanglement exhibits vari- 
ous new features. Notably, there are two different equiv- 
alence classes of genuine three-qubit entanglement 
the representatives being any one of the two maximally 
entangled Greenberger-Horne-Zeilinger (GHZ) states [2j 
|GHZ ± ) = (|TTt) ± |U4))/\/2 on the one hand, and the 
W state i |W) = (|TTI> + ITIT) + IITT»/v / 3 on the other. 
The ability to realize both representatives in real phys- 
ical systems is thus of high importance in the study of 
genuine tripartite entanglement. Particularly interesting 
is the GHZ state, as it represents the strongest quantum 
correlations possible in a system of three qubits. Further- 
more, it is equivalent to the three-qubit cluster state used 
in one-way quantum computation It is favorable to 
obtain the GHZ and W states as the eigenstate of a suit- 
able system, rather than by engineering them using quan- 
tum gates. In this Letter, we demonstrate the possibility 
of obtaining approximate GHZ and exact W states as the 
ground state (g.s.) of three spin-qubits in a ring geometry 
coupled via an anisotropic Hciscnbcrg interaction. The 
use of quantum gates is therefore not required. Rather, 
the desired states are achieved merely by cooling down to 
sufficiently low temperatures. We state all our results in 
terms of the exchange coupling strengths in order to keep 
our proposal open to a broad set of possible implemen- 
tations of the qubits. We remark that, while Heisenberg 
models have been studied frequently in the context of 
entanglement Q (also with respect to entangled eigen- 
states 0), this is the first time that highly entangled 
states are reported as the non-degenerate g.s. of three 
exchange-coupled qubits. Our study inevitably involves 
the issue of quantifying entanglement 
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At finite temperatures, the mixing of the g.s. with excited 
states forces us to evaluate a mixed-state entanglement 
measure (EM) in order to study the entanglement in the 
system meaningfully. Computationally, this is a rather 



formidable task. We generalize a numerical scheme that 
has originally been developed to compute the entangle- 
ment of formation (EOF) 0, Il3j |. Our scheme can be 
used to evaluate any mixed-state EM defined as a so- 
called convex roof hj. 

Model. — We assume that three spins Sj, with S = 1/2, 
are located at the corners i = 1, 2, 3 of an equilateral tri- 
angle lying in the .ry-plane. Their interaction is described 
by the anisotropic Heisenberg Hamiltonian 
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H = — J X y^^(Sf Sf +1 + SfSf +1 ) — J z ^ SiSf+i + Hz, 



i=l 



i=l 



(1) 

where S4 = Si. Here, J xy and J z are the in- and out- 
of-planc exchange coupling constants, respectively, and 
Hz = Y^i=i k» ' denotes the Zeeman coupling of the 
spins Sj to the externally applied magnetic fields bj at the 
sites i [l5| • We now seek a configuration of b; 's yielding 
a highly entangled GHZ- or W-type ground state. Finite- 
temperature effects will then be studied in a second step. 
Ground-state properties. — We first consider isotropic 



exchange couplings, i.e., J xy 



J z = J. For b; 



0, we 



naturally find two fourfold-degenerate cigenspaces due 
to the high symmetry of the system. For J > 0, i.e., 
ferromagnetic coupling, the ground-state quadruplet is 
spanned by the two GHZ states |GHZ ± ), the W and 
the spin-flipped W state. Appropriately chosen magnetic 
fields allow one, however, to split off an approximate GHZ 
state from this degenerate eigenspace. To identify the op- 
timal field geometry, we first observe that the two states 
|GHZ ± ) have the form of a tunnel doublet. If we thus find 
a set of bj's, which, in the classical spin system, results 
in precisely two degenerate minima for the configurations 
t and HI with an energy barrier in between, quantum 
tunneling will yield the desired states. In order to sin- 
gle out exactly the two directions perpendicular to the 
xy-plane, the magnetic fields must be in-plane, be of the 
same strength, and sum to zero. This immediately im- 
plies that successive directions of the fields must differ by 
an angle of 27r/3 from each other. We choose the fields 
to point radially outwards, although any other config- 
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uration possessing the required symmetry is equivalent. 
However, this setup is experimentally most feasible, e.g., 
by placing a bar magnet below the center of the sample 
(in the case of a solid state implementation) . In order to 
favor parallel spin configurations we consider the regime 
where J S> b, b = |b,| being the Zeeman energy. We 
may thus assume that for given mean spherical angles 
■d (zenith) and (p (azimuth), the orientation of each spin 
will deviate from these values only by a small amount. 
Expanding the classical energy E c ($, (p) corresponding to 
Eq. ([1]) to second order in these deviations and minimiz- 
ing with respect to them under the constraint that they 
separately sum to zero yields: 
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sin (3<^)sin 3 tf. (2) 



This expression is minimal for $ = and $ = 7r, rep- 
resenting the desired configurations. The paths in i9 
with lowest barrier height connecting these two minima 
are found for values of (p = + 2irn/3 mod 2tt, 

n = 0, 1, 2, reflecting the rotational symmetry of the sys- 
tem. The corresponding barrier height is approximately 
given by [(6/ J) 2 - (&/J) 3 /6]/4 [lj]. 

Next we return to the quantum system. The Hamil- 
tonian (fT]) with isotropic exchange coupling J and ra- 
dial magnetic field can be diagonalized exactly. Expand- 
ing for bj J <C 1, the overlap probabilities of the exact 
ground state |0) with |GHZ + )i. u . and the exact first ex- 
cited state |1) with |GHZ~) l u ., respectively, are identi- 
cal to second order and are given by |i. u .(GHZ + |0)| 2 = 
| Lu .(GHZ-|l)| 2 w 1 - |(6/ J) 2 ('l.u.' indicates that 
the states are equivalent to GHZ states via local uni- 
tary transformations). The associated energy splitting is 
given by AE ,i ~ 2(6/J) 3 /3 (sec inset of Fig. EE). This 
confirms the above semiclassical considerations in terms 
of tunnel doublets. Moreover, we see that the g.s. can 
only approximate a GHZ state although this approxima- 
tion will turn out to be very good even at finite tem- 
peratures where mixing with excited states additionally 
decreases the entanglement. Before discussing this in 
greater detail, we study the ground-state of the general 
anisotropic case with J xy ^ J z in the Hamiltonian (TIJ. 

When J xy 7^ J z it is possible to generate highly en- 
tangled states by applying a spatially uniform magnetic 
field either perpendicular to or in the xy-plane. Indeed, a 
field along the z-axis, i.e., hi = be z , i = 1,2, 3, with b > 
yields an exact W state as g.s. if J xy > and b < J xy — J z 
(note that this implies the condition J xy > J z ). The 
optimal Zeeman energy 6 opt leading to the highest en- 
ergy splitting AE opt between the g.s. and the first ex- 
cited state is given by 6 pt = {Jxy — Jz)/2- This yields 
A£ opt = 3J xy /2 if J z < -2J xy and A£ opt = (J xy -J z )/2 
otherwise. The W state is thus best realized by choosing 
b = b opt together with a temperature sufficiently small 
compared to A_E opt . In order to obtain a GHZ state, 
one has to apply an in-plane magnetic field b; = be x . In 



this case we find for J z > 0, — 2 J z < J xy < J z a situation 
similar to the one in the case of isotropic coupling and 
radial magnetic field: The g.s. converges to a GHZ state 
for vanishing field but also the energy difference to the 
first excited state goes to zero in this limit. 

Entanglement measure. — Below, we will quantita- 
tively study the effects of finite temperature T > on 
the amount of entanglement present in the system. For 
this purpose, we evaluate a suitable mixed-state EM of 
the canonical density matrix p of the system. The three- 
tangle, or simply tangle r p (originally called residual en- 
tanglement), is an EM for pure states £ Tti(£)Tt2 (8)0-3 
of three qubits. It reads [11 1 



T v {\i>)) =4detTr 2 , 3 Pp-C 2 (Tr 3 p p )-C 2 (Tr 2j0p ), (3) 

where p p = \ip) (ip\ , Tr^ denotes the partial trace over sub- 
system i, and C is the two-qubit concurrence The 
tangle takes values between and 1 and is maximal for 
GHZ states. It is also known that t p is an entanglement 
monotone The generalization of pure-state mono- 
tones to mixed states is given by the so-called convex 
roof 0, [3, Ell • Accordingly, the mixed-state tangle r is 
defined as 



r( P ) 



inf 

{p«,|Vi>}es(p) 



(4) 



Here, D(p) denotes the set of all pure-state decomposi- 
tions {pi, |V>t)}iLi of p, with 2i=lft|V'i}('0i| = Pi > °) 
1 and K > R = rank p. The above definition 



of t ensures that r{p) = T p (\ip)) if p = \ij;){tp\, and that 
r itself is an entanglement monotone Q . 

Numerical evaluation. — In order to tackle the opti- 
mization problem in Eq. (j4|) numerically, the set of all 
pure-state decompositions 'D(p) needs to be given in an 
explicitly parameterized form. It is known 20, 21 1 that 



every pure state decomposition {pi, \ipi)}fLi °f P 1S rc " 
lated to a complex K x R matrix U satisfying the unitary 
constraint Wll = lowj, i.e., a matrix having orthonor- 
mal column vectors [22|. In fact, the set of all such ma- 
trices, the so-called Stiefel manifold St(K,R), provides 
a complete parametrization of all pure-state decompo- 
sitions {pi, \ipi)}jL 1 £ D(p) of p with fixed cardinality 
K. The minimization problem in Eq. ([4]) can thus be 
rewritten as 



f(p)= min inf h(U.p). 
Jxr ' K>RuaSt(K,R) y ,Hh 



(5) 



where in our case h is the sum over the weighted pure- 
state tangles with probabilities and state vectors ob- 
tained from p via the matrix U. Problems of this kind 
are considered to be extremely difficult to solve in gen- 
eral Q. We have performed the minimization over the 
Stiefel manifold numerically using the method described 
below. We have found that the thereby obtained values 
converge quickly as K is increased, and have thus fixed 
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K = R + 4 throughout all of our calculations, yielding 
an accuracy which is by far sufficient for our purpose 
(note that decompositions with smaller cardinality are 
contained as well). The numerical method we used is 
a generalization of the conjugate gradient algorithm pre- 
sented in Ref. [T3| ■ It is however only suited for searching 
over the unitary manifold St(K, K). At the cost of over- 
parameterizing the search space, we have to minimize 
over K x K matrices using only the first R columns. The 
iterative algorithm builds conjugate search directions X 
(skew-hcrmitian K x K matrices) from the gradient G 
at the current iteration point U and the previous search 
direction using a modified Polak-Ribiere update formula. 
A line search along the geodesic g(t) = U exp(tX) going 
through U in direction X is performed in every step. In 
Ref. [l3j, an analytical expression for the gradient G is 
given in the case where / is the EOF. The algorithm is 
however also applicable to a generic convex-roof EM / 
of the form ((5]). We find the matrix elements Gjk of the 
general gradient G to be 

G jk = (A jk - A kj )/2 + i(S jk + S kj )/2, (6) 

where 

^=E(^k Rc ^ + ailk Im ^)' (7) 

i—l v / 

^ = EfTrV Rc ^-^rV Im ^V ( g ) 

^ \dlmU ik dReUi k J J 

The derivatives of h with respect to the real and imag- 
inary parts of Uik, Ret/jfc and ImC/^., respectively, are 
taken at U and can be evaluated numerically using fi- 
nite differences. We have tested our implementation 
by comparing our numerical results to known analyti- 
cal results. The maximal encountered absolute error was 
smaller than 10~ 13 for the EOF of isotropic 2x2 states 
1CT 12 for 3 x 3 states and 1(T 10 for the tangle of a 
GHZ/W mixture This suggests that, although our 
method can only provide an upper bound, this bound is 
very tight. It was shown only recently that also a (typi- 
cally tight) lower bound on any entanglement monotone 
can be estimated using entanglement witnesses 0, [HI . 
This is an interesting subject which is left for future re- 
search. 

Finite temperature. — We return to the study of the 
three qubits described by the Hamiltonian {T]). Using the 
generalized conjugate gradient algorithm, we are able to 
investigate the entanglement as a function of the temper- 
ature T, the magnetic field strength b and the exchange 
couplings J xy and J z by calculating the mixed-state tan- 
gle r(p), where p = exp(-H/k B T)/ Tr e^q>{- H / k B T) is 
the canonical density matrix of the system. To our knowl- 
edge, this is the first time that r(p) has been evaluated 
for states arising from a physical model. Our main goal 
now is to maximize the entanglement as a function of 
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FIG. 1: (color online) The tangle r of the system with 
isotropic positive (ferromagnetic) coupling J and radial mag- 
netic field as a function of b/J for different temperatures 
T = l(r 4 J/fc B (dashed line), lCr 3 J/fc B (dash-dotted line), 
10 _2 J/fcB (dash-dot-dotted line) and 5 x 10 _2 J/fcB (dotted 
line). Note the twofold influence of the temperature on r: 
Although higher temperatures reduce the maximally achiev- 
able entanglement, a stabilizing effect is observed as well. A 
maximum in the tangle is more robust against fluctuations 
in b at higher temperatures due to the less rapid drop-off 
of t as b/J is reduced. Conversely, r of the approximate 
GHZ g.s. |GS) (T = 0, red line) shows a discontinuity at 
6 = 0, where r(p) = 0. For b > 0, we find the simple al- 
gebraic expression r p (|GS}) = (3 — 8b/J)/C + 2/y/C, where 
C = 9 + 4b{Ab/J - 3) /J. Inset: Energy splitting AE ,i of the 
ground-state doublet as a function of b/J. 



b = |bj|, i.e., the Zeeman energy. For this purpose we 
consider only GHZ states in the following, since our W 
ground states are b- independent (see above). In the sys- 
tem with isotropic exchange coupling J > and radial 
magnetic field, the tangle r tends to zero for b/J^0 
due to the vanishing energy splitting Ai^o.i (see Fig.[T|). 
We remark that this behavior is discontinuous at T = 0, 
where r(p) -* 1 for b/J -> 0, but r(p) = at b = 0. 
With larger b/J, the g.s. contributes dominantly to p 
but simultaneously deviates increasingly from a GHZ 
state. The entanglement in the system is therefore re- 
duced (cf. solid line in Fig. [I}. For a given temperature, 
the maximal tangle r max is therefore obtained at a fi- 
nite optimal value (6/J) pt of the scaled magnetic field 
strength as a trade-off between having a highly entangled 
g.s. and separating the latter from excited states in order 
to avoid the negative effects of mixing. For low tempera- 
tures T < 10 -2 J/fce, we numerically find the power laws 
(6/J) op t oc (k B T/J) a and 1 - r max c< {k B T/jf with the 
exponents a w 0.30 and (3 « 0.63. Specifically, we ob- 
tain t( P ) = 0.98 (0.92) for T = Ur A J/k B (10~ 3 J/k B ) 
and b = 0.11J (0.21J). Apart from the effect of reducing 
Tmax, finite temperatures also possess the advantageous 
feature of broadening the discontinuity of r at T = and 
b = which makes r max more stable against fluctuations 
of b around 6 op t (see Fig. [1]). 

We finally come back to the general anisotropic 
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FIG. 2: (color online) Top: Maximally achievable tangle r max 
in the anisotropic system (GHZ g.s.) with homogeneous in- 
plane magnetic field and J z > as a function of temperature 
for six anisotropy ratios J xy /J z (see legend). The curves end 
at r(p) = 10 -5 . Bottom: The corresponding optimal values 
(6/J 2 ) op t of the scaled magnetic field strength b/J z . 



model ([T]) with J xy ^= J z subject to a homogeneous in- 
plane magnetic field. In Fig. [5] we show the maximally 
achievable tangle r max (optimized with respect to b/J z ) 
as a function of temperature for various anisotropy ratios 
Jxy/ Jz (where, as before, J z > 0). Since we are inter- 
ested in high values of T max , an arbitrary but low cutoff 
was introduced in the calculation at r{p) = 10~ 5 . The 
lower panel of Fig. [2] depicts the corresponding optimal 
field values (b/J z ) op t- At low temperatures T, a power- 
law dependence of (b/ J z ) opt on T is observed, similar to 
the above isotropic case. Note that a higher amount of 
entanglement can be realized in systems with stronger 
anisotropics. E.g., for Ising coupling (J xy /J z = 0) we 
find t( P ) = 0.98 (0.89) for T = lQ- A J/k B (10~ 3 J/fc B ) 
and b = 0.080J Z (0.16J 2 ). At T = W^J/kB but 
with Jxy/Jz = 0.9, still a very good value r(p) = 0.90 
is achieved for 6 = 0.016J Z . We remark that still 
higher tangles arc obtained for negative (antiferromag- 
netic) J xy > —2J Z . In this case, the maximal tangle as 
a function of T decays even more slowly than the curves 
displayed in the top panel of Fig. O 

Possible implementations of the qubits include GaAs 
and In As quantum dots, In As nanowires or single- wall 
carbon nanotubes. Assuming a typical value of \J\ ~ 
lmeV HI, we obtain r ft 0.9 at T ft 10 mK and 
B ft 2T (assuming a g-factor of \g\ = 2). Ferromagnetic 
coupling is achieved by operating the dots with more than 
one electron per dot. 
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